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Abstract 

We review some recent results beyond the now established theory of magnetiza- 
tion switching of a nanoparticle within the single-spin approximation. The first 
extension is that of the Stoner-Wohlfarth model for magnetization static switching 
under applied magnetic field including the effect of temperature at long-time scales. 
The second concerns a generalization of the Neel-Brown model for thermoactivated 
dynamic magnetization switching to include the effect of exchange interaction in 
the framework of Langer's theory in the intermediate-to-high damping limit. We 
finally argue why the single-spin approximation is not appropriate for very small 
nanoparticles. 
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1 Introduction 



Magnetic nanoparticles are very interesting and challenging systems to both 
fundamental research in physics and technological applications. From the point 
of view of fundamental physics, nanoparticles are interesting because they offer 
a rich laboratory for the study of dynamic phenomena as they show super- 
paramagnetism at high temperature and exponentially slow relaxation at low 
temperature. On the other hand, due to their high coercivity, their magnetiza- 
tion enjoys long-range stability, a very important property for the information- 
storage technology. In principle, particles of very small size may be used to 
increase the storage density. However, this is hampered by many difficulties 
related with the finite size of the particles. In particular, surface and thermal 
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effects become dominant in such particles, and drastically affect the magne- 
tization relaxation time, and thereby the stability of the information stored 
in recording media. Therefore, the aim of any study of such systems, be it 
experimental or theoretical, is to understand the dynamic properties of small 
nanoparticles taking account of spatial inhomogeneities due to, e.g., surface 
disorder. For this purpose one has to tackle the difficult problem of including 
surface effects in the calculation of the relaxation time of the particle magneti- 
zation. However, this requires a microscopic approach to account for the local 
environment inside the particle, and thus include microscopic interactions such 
as spin-spin exchange, together with magnetocrystalline anisotropy and Zee- 
man interaction. Unfortunately, this leads to a rather difficult task owing to 
the large number of degrees of freedom which hinders any attempt to analyze 
the energyscape. For this reason, inter alia, calculations of the reversal time 
of the magnetization of fine single-domain ferromagnetic particles initiated by 
Neel [1], and set firmly in the context of the theory of stochastic processes by 
Brown [2] , [3] , have invariably proceeded by ignoring all kind of interactions. 
Thus, the only terms which are taken into account in these calculations are 
the internal magneto-crystalline anisotropy of the particle, the random field 
due to thermal fluctuations, and the Zeeman term. This assumes that the 
particle's atomic moments (or spins) switch in a coherent manner, so that the 
magnetic state of the particle can be described by a giant magnetic moment, 
this is the one-spin approximation, and the particle is then dealt with as a 
one-spin system. 

It is well known that the magnetization of a nanoparticle can overcome the 
energy barrier, due to magneto- crystalline anisotropy, and thus reverse its di- 
rection, at least in two ways ^ : either under applied magnetic field which 
suppresses the barrier, or by thermal effects which produce statistical fiuctua- 
tions. Within the framework of the one-spin approximation, the magnetization 
switching under applied magnetic field, at zero temperature, is well described 
by the Stoner-Wohlfarth model [5]. This model has been confirmed by ex- 
periments on single cobalt particles by Wernsdorfer et al. [see Ref. [6] for a 
review]. At finite temperature, but at long-time scales or quasi-equilibrium, 
this switching occurs according to two regimes. At very low temperature, this 
is due to the coherent rotation of all spins, as in the Stoner-Wohlfarth model, 
whereas at higher temperature, the magnetization switches by changing its 
magnitude ^ . This results in a shrinking of the Stoner-Wohlfarth astroid as 
predicted by the modified Landau theory [7], and confirmed by experiments 
[6]. At short-time scales, and within the one-spin approximation, crossing of 



^ It has quite recently been shown [4], experimentally and theoretically, that effi- 
cient magnetization switching can be triggered by transverse field pulses of a dura- 
tion that is half the precession period. 

^ Rigorously, the switching of magnetization through a change of its magnitude 
cannot, however, be explained within the single-spin approximation [see section 4]. 
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the energy barrier activated by thermal energy is described by the Neel-Brown 
model [1], [2], [3] and its extensions [see section 2.3], again in agreement with 
experiments on individual cobalt particles [6]. 

In both Stoner-Wohlfarth and Neel-Brown models, the energy of a fine single- 
domain ferromagnetic nanoparticle only contains the internal magneto-crystalline 
anisotropy and the Zeeman term. The thermoactivated dynamics of the par- 
ticle, e.g., within the stochastic Langevin approach, are described by adding 
a stochastic field in the equations of motion of the magnetization, in order to 
account for thermal effects exerted on the particle by the surrounding bath. 
This stochastic field is responsible for the statistical fiuctuations of the mag- 
netization direction. 

In this paper we briefiy describe these two models. Then, we present two ex- 
tensions thereof: the first is of the Stoner-Wohlfarth model to include thermal 
effects on the magnetization switching at long-time scales, or equivalently, 
at quasi-equilibrium. The second extension is that of the Neel-Brown model 
which takes account of exchange interaction and investigates its effect on the 
relaxation time of a pair of magnetic moments, within Langer's approach. It 
will be shown that the Neel-Brown result for the relaxation time of a rigid 
magnetic moment is only recovered for very large exchange coupling, while 
for weak coupling there appear new interesting features, giving a foretaste of 
the intricacies and subtleties that should arise in a particle when treated as a 
multi-spin system. 

The paper is organized as follows: In section 2 we deal with a nanoparticle in 
the single-spin approximation. For such a single-spin system, we briefiy discuss 
the Stoner-Wohlfarth model describing the magnetization switching under the 
applied magnetic field and then consider the effect of temperature within the 
framework of Landau's theory. Next, the Neel-Brown model for thermoacti- 
vated switching is reviewed. We comment on the applicability of this model, 
and its further extensions to non-axially symmetric potentials, to the calcula- 
tion of the blocking temperature of an assembly of nanoparticles. In section 3, 
we abandon the single-spin approximation and deal with the simplest, though 
non trivial, problem of two exchange-coupled spins and study the effect of the 
latter on the relaxation time within Langer's approach. The results are then 
compared with those of the Neel-Brown model. The last section motivates and 
discusses the generalization of the latter work to a multi-spin particle. 



2 Magnetization switching of a nanoparticle as a single-spin system 

Before we proceed, we would like to stress the role of anisotropy. Exchange 
energy is completely isotropic, that is it does not depend on the direction in 
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space in which the crystal is magnetized. One often chooses the reference z 
direction as the one along which the component of the magnetization is calcu- 
lated. However, this direction does not really have any meaning in the limit of 
zero applied field. In fact, due to global rotation the net magnetization van- 
ishes in zero applied field. This not only contradicts experiments but it is also 
in confiict with everyday experience that, for instance, the particles in an au- 
dio or video tape remain magnetized, and do not lose the recorded information 
upon switching off the writing field. This is so because real magnetic materi- 
als are not isotropic, and that the reference z direction in question is actually 
defined by anisotropy. There are several types of anisotropy, but the most 
common is the magneto-crystalline anisotropy, caused by the spin-orbit inter- 
action. Indeed, the electron orbits are linked to the crystallographic structure, 
and by their interaction with the spins they make the latter prefer to align 
along well-defined crystallographic axes. There are therefore some directions 
in space in which it is energytically more favorable to magnetize a given crystal 
than in other directions, the difference being given by the direction-dependent 
anisotropy energy. As for the exchange interaction, quantitative estimates of 
the spin-orbit interaction from basic principles are also possible but the ac- 
curacy is not as good. Therefore, the anisotropy term is always written as a 
phenomenological expression, which is actually the first term of a power series 
expansion that take into account the crystal symmetry, and the corresponding 
coefficient is usually taken from experiment. The magneto-crystalline energy 
is usually 2 — 4 orders of magnitude smaller than exchange energy, but the 
direction of the magnetization is determined only by anisotropy. On the other 
hand, exchange coupling tends to align all the spins inside the particle parallel 
to each other, while anisotropy tends to align them along a certain crystallo- 
graphic direction. A compromise between exchange and anisotropy is obtained 
by aligning all spins parallel to each other and to the anisotropy direction. 



Other features of a magnetic material, namely its subdivision into domains, 
is then caused but yet another energy term, called the magnetostatic self- 
energy (demagnetization or shape energy). This is the surface term that stems 
from the dipolar interactions. It is irrelevant for spherical samples. However, 
for many materials, slight deviation from spherical shape renders the shape 
anisotropy very important especially for the spatial dependence of the magne- 
tization. On the other hand, as shown by Brown and Morrish in [8], a single- 
domain particle with an arbitrary shape is equivalent to a suitably chosen 
general ellipsoid of the same volume, as far as the total energy is concerned. 
In all the sequel, we shall assume that the particle under study is spherical 
and single-domain. 
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2. 1 Magnetization switching under magnetic field: Stoner- Wohlfarth model 



In small magnetic particles exchange interaction may be strong enough to 
hold all spins tightly parallel to each other, and prevents spatial dependence 
of the particle's magnetization. Wc have shown [7], [9] that this no longer holds 
when free-boundary effects and/or surface anisotropy arc taken into consider- 
ation. Ignoring surface effects, or in other words spatial inhomogeneities, the 
exchange energy is a constant and plays no role in the energy minimization. In 
the latter only the anisotropy energy and the interaction with the applied field 
are relevant. Accordingly, let us consider the simplest case of a magnetic mo- 
ment m with a uniaxial anisotropy axis e along the z direction of the applied 
field. The corresponding Hamiltonian reads 

Kv 

-(m.e)2-m.H, (1) 



where K is the magneto-crystalline anisotropy constant and v is the volume 
of the particle. Upon writing m = ms, H = Heh, and introducing the dimen- 
sionless anisotropy and field parameters 



Kv 



h 



Hm 



(2) 



where ks is the Boltzmann constant, the Hamiltonian (1), divided by thermal 
energy, becomes 



-pn^a {s.ef + 2h{s.eh) 



(3) 



Denoting by 9 the angle between the magnetization direction s and the field 
direction Sh, we write (3) as 



pn 



cos^e + 2/tcos6' 



(4) 



For simplicity, we only consider the case of easy axis, i.e., a < 0. For o" > 
(easy plane) the results are the same only that the energy maximum and 
minimum are interchanged. The extrema of and their nature are obtained by 
setting to zero its first derivative with respect to 9, and evaluating the second 
derivative at these extrema. The results are presented in Table 1. Thus, for 
\h\ < 1 the energy has minima at ^ = and 6 = n, separated by a maximum at 
9„i = arccos(— /i). For \h\ > 1 the upper (also the shallower) energy minimum 
{9 = 71 for h > 0) turns into a maximum as it merges with the intermediate 
maximum at 9^, which disappears (see Fig. 1). From the values of the energy 
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Table 1 



Field 


Minima 


Maxima 


\h\<l 


61 = 0,7r 


9 = arccos(— /i) 


h>l 


^ = 


e = Tr 


h<-l 


e = Tr 


6 = 



PH{e) 






/? = / / 




^ h=0.25y / 




\h=0.75/ 





0.00 0.25 0.50 0.75 g/^ 1.00 



Fig. 1. Magnetic energy in the case of longitudinal field for some values of the 
reduced field h. Upon increasing the field the number of potential wells changes 
from two to one. 




-1.0 -0.5 0.0 0.5 , 1.0 -1.0 -0.5 0.0 0.5 , 1.0 

h h 



Fig. 2. Left: (numerical) hysteresis loops for different values of ip increasing inwards: 
tp = 0, 60°, 85°, 90°, for a 3^ particle with uniaxial anisotropy and strong exchange 
coupling. Right: (numerical in squares and analytical in full line) SW astroid for 
the same particle. 
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at ^ = 0, TT and, when it exists, at the intermediate maximum 6m, one obtains 
the energy-barrier heights {\h\ < 1) 



AH^P [Hiem) - n{0, tt)] = (7(1 ± hf. 



(5) 



The analysis above helps solve for the hysteresis curve, an important part of 
magnetism, for small ferromagnetic particles. Such calculations are known as 
the Stoner- Wohlfarth model [5] . In their original study, Stoner and Wohlfarth 
considered shape anisotropy instead of uniaxial anisotropy. However, their 
model is usually applied to the latter. In addition, they also considered the 
more general situation of a field applied at an arbitrary angle ip to the easy 
axis. Prom the above analysis we see that there is a unique minimum for > 1 
while there are two minima for \h\ < 1. This is due to the multivaluedncss of 
the trigonometric functions entering the Hamiltonian (4) and its derivatives 
with respect to 9. In order to obtain a unique solution, one has to specify 
and follow the history of the field h for each angle ip. One usually starts at 
saturation and increases (or decreases) the field across zero until it reaches 
the value at which the energy barrier disappears. This field marks the limit of 
metastability and is called the critical field. In Stoner- Wohlfarth model, this 
is given by = 1 [for ip = 0]. For arbitrary ip the critical field is given by 



One can also define what is called the switching field, that is the field at which 
the magnetization changes sign. 

In Fig. 2 we present (on the left) the hysteresis loop at different angles ip, and 
on the right (in full line) the angular dependence of the switching field, the so- 
called Stoner- Wohlfarth astroid, as obtained from the Stoner- Wohlfarth model 
[see Eq. (15) with a = 0]. In fact, the hysteresis loops in Fig. 2 and the astroid 
in squares, were obtained from the numerical solution of the Landau-Lifshitz 
equation [9]. 

Stoner- Wohlfarth astroids for cubic anisotropy have been obtained by Thiaville 
[10] using a geometrical approach. 

2.2 Thermal effects on the Stoner- Wohlfarth model: Landau theory 

Before dealing with the effect of thermoactivation on the magnetization switch- 
ing between the different energy minima, and discuss the calculation of the 
corresponding relaxation time, we now consider the effect of temperature on 
the Stoner- Wohlfarth model [7]. More precisely, we consider magnetization 



he = 



1 




3/2 ■ 
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switching at a long-time scale, i.e., at quasi-equilibrium. In [7] wc derived a 
free energy for weakly anisotropic fcrromagnets which is valid in the whole 
range of temperature and interpolates between the micromagnetic energy at 
zero temperature and the Landau free energy near the Curie point Tc. This free 
energy takes into account the change of the magnetization magnitude due to 
thermal effects, in particular, in the inhomogeneous states. As an illustration, 
we studied the thermal effect on the Stoner-Wohlfarth astroid and hysteresis 
loop of a ferromagnetic nanoparticle assuming that it is in a single-domain 
state. Within this model, the saddle point of the particle's free energy, as 
well as the metastability boundary, are due to the change in the magnetiza- 
tion magnitude sufficiently close to Tc, as opposed to the usual homogeneous 
rotation process at lower temperatures. 

For weakly anisotropic magnets, where the anisotropy energy is much weaker 
than the homogeneous exchange energy, the magnetization magnitude M is 
either small or only slightly deviates from Me, the equilibrium magnetization 
at zero field. Thus the condition for the magnetization minimizing the free en- 
ergy derived in [7] , \M — Mg \ -C Mg , where Mg is the saturation magnetization 
at T = 0, is satisfied in the whole temperature range. In [7] we presented the 
derivation of this free energy and illustrated some of its features in the case 
of a single-domain magnetic particle with a uniaxial anisotropy. 

For single-domain magnetic particles in a homogeneous state, the gradient 
(exchange) terms in the free energy can be dropped and the free energy can 
be written in the form 



F = + {VM'Jx±)f 

/ = -n-h+^(n^ + nJ) + ^(n2-lf, (6) 

where V is the particle's volume and / the reduced free energy in terms of the 
reduced variables 

n = M/Me, h = Hx±/Me, a = 2x\\/x±- (7) 



where x\\/ X± the longitudinal and transverse susceptibilities. We see 
that the parameter a here controls the rigidity of the magnetization vector; it 
goes to zero in the zero-temperature limit (with fixed magnetization modulus) 
and diverges at Tc. 

For fields h inside the Stoner-Wohlfarth astroid [see Fig. 2 (right)], which will 

be generalized here to nonzero temperatures, / has two minima separated by a 
barrier. Owing to the axial symmetry, one can set Uy = for the investigation 
of the free energy scape. The minima, saddle points, and the maximum can 
be found from the equations df /dux = df /duz = 0, or, explicitly 
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Fig. 3. The free energy of a ferromagnetic particle with uniaxial anisotropy [/ in 
Eq. (6) in zero field] for a = 2x||/x± = 0-5 (left) and a = 2 (right) corresponding to 
lower and higher temperatures, respectively. 

77,^(71^ - 1) = ah:, 

nx{n^ — 1 + a) = ahx- (8) 
Solving for Ux one obtains a 5th-order equation for 

hlnl^ {h, + n,y{aK + n,-nl). (9) 



In zero field, the characteristic points of the energyscape can be simply found 
from Eqs. (8). One of these points is Ux = Uz = 0, which is a local maximum 
for a < 1 and a saddle point for a > 1. The minima are given by Ux — 0, 
Uz = ±1. The saddle points correspond to Uz — 0, while from the second of 
Eqs. (8) one finds 



nx = l ' ~ (10) 




In fact, due to the axial symmetry, for a < 1 one has a saddle circle n'^ + Tiy = 
1 — a, rather than two saddle points. The free-energy barrier following from 
this solution is given by 

' (2 - a)/4, a < 1 

A/ = /sad - /min = r " (H) 

l/(4a), a > 1. 



The free-energy landscape in zero field is shown in Fig. 3. At nonzero tem- 
peratures a > 0, the magnitude of the magnetization at the saddle is smaller 
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than unity since it is directed perpendicular to the easy axis, and for this 
orientation the "equihbrium" magnetization is smaller than in the direction 
along the z axis. For a > 1, the two saddle points, or rather the saddle circle, 
degenerate into a single saddle point at n^, = = 0, and the local maximum 
there disappears. That is, for the magnetization to overcome the barrier, it 
is easier to change its magnitude than its direction. This is a phenomenon of 
the same kind as the phase transition in ferromagnets between the Ising-like 
domain walls in the vicinity of T^. (the magnetization changes its magnitude 
and is everywhere directed along the z axis) and the Bloch walls at lower 
temperatures [11], [12]. 

The Stoner-Wohlfarth curve separates the regions where there are two minima 
and one minimum of the free energy. On this curve the metastable minimum 
merges with the saddle point and loses its local stability. The corresponding 
condition is 

d'^fjdnl X d'^fldnl - {d^f/dn^dn.f = 0, (12) 

which upon using Eq. (9) leads to the quartic equation for 

h^[{2 + a)n^ + 3ah^] + 2nl ^ 0. (13) 

Before considering the general case, let us analyze the limiting cases a <^ 1 
and a ^ 1. 

At low temperatures, i.e., a <^ 1, the magnetization only slightly deviates 
from its equilibrium value, and to first order in a, one obtains 

n^^\- an\ or + (1 + a)nl ^ 1. (14) 

From Eq. (13) and the analogous equation for n^; one derives the equation for 
the Stoner-Wohlfarth astroid 

/i^/=^ + [(l + a/2)/i,]2/3^1, a«l. (15) 

One can see that, in comparison with the standard zero-temperature Stoner- 
Wohlfarth astroid, i.e. at a = 0, is rescaled. The critical field in the z 
direction decreases because of the field dependence of the magnetization mag- 
nitude at nonzero temperatures. 

Similarly, in the case a ^ 1, i.e. near Tc, Eq. (12) leads to the equation for 
Stoner-Wohlfarth curve 

nl + 3nl^l. (16) 
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Fig. 4. Left: Dependence hz{a) at hx = on the Stoner-Wohlfarth curve. Right: 
The Stoner-Wohlfarth curves at different temperatures, a = 2x\\/x± = = 0); 
0.3, 0.5, 2/3, 0.835, 1, 2. 

and upon using Eq. (8) we infer 

^ -{ah,/2f'\ ^ h,. (17) 



Making use of this result in Eq. (16), one obtains another limiting case of the 
Stoner-Wohlfarth curve 

hi + 3(a/i,/2)2/3 ^1, a > 1. (18) 



In this case, the critical field (i.e., the field on the Stoner-Wohlfarth curve) in 
the z direction is strongly reduced, and there is no singularity in the depen- 
dence hczihx) ai hx — 0. 

The quahtatively different character of the Stoner-Wohlfarth curves in these 
two cases is due to the different mechanisms pertaining to the loss of the local 
stability for the field applied along the z axis. For hx — the mixed derivative 
{d"^ f / duxdriz) vanishes and Eq. (12) factorizes, 

{a-l + nl){-l + 3nl) =0. (19) 



Vanishing of the first factor in this equation corresponds to the loss of stability 
with respect to the rotation of the magnetization, d'^f/dnl = 0, and vanishing 
of the second factor, d'^f/dnl = 0, imphes the loss of stabihty with respect to 
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the change of the magnetization magnitude. Using Ux — 0, with the help of 
the first equahty in Eq. (9) one obtains in the two cases 



Note that the transition between the two regimes occurs here at a different 
value of a than in Eq. (10). The dependence hz{a) at /i^ = is shown in Fig. 4 



In the general case, it is easier to find the Stoner-Wohlfarth curve numerically 
from Eqs. (8) and (13). The results in the whole range of a are shown in Fig. 4 
(right). 

We have shown that, for single-domain magnetic particles with uniaxial anisotropy, 
thermal effects qualitatively change the free-energy landscape at sufficiently 
high temperatures, so that the passage from one free-energy minimum to the 
other is realized by the uniform change of the magnetization magnitude rather 
than the uniform rotation. This also qualitatively changes the character of 
the Stoner-Wohlfarth curve and hysteresis loops [7]. These effects cannot be 
observed with standard methods, however, because keeping the height of the 
free-energy barrier much larger than thermal energy requires so large particle 
sizes that the single-domain criterion is no longer satisfied [see discussion in 
Sect. 6 of Ref.[7]]. For the uniform states, the theory is valid at low temper- 
atures, but then the thermal effects considered here are small corrections to 
the zero-temperature results. 

Finally, we would like to mention that quite recently, experimental results have 
been obtained by Jamet et al. [13] on 3 nm cobalt nanoparticles which clearly 
show the disappearance of the singularity near if j; = at a temperature circa 
8 K (the blocking temperature being 14 K). The height of the experimental 
astroid decreases nearly as its width with increasing temperature, but it does 
not become flat as predicted by our calculations. This failure of our theory, 
in addition to the fact that this effect is predicted at much higher tempera- 
tures, is not surprising considering the fact that T <^ Tc. Nevertheless, the 
disappearance of the singularity is definitely predicted by our theory. 

2.3 Magnetization switching as a thermal effect: N eel-Brown model 

We have seen that the Stoner-Wohlfarth model accounts for the hysteretic 
rotation of the particle's magnetization over the potential barrier under the 
infiuence of a field applied in an arbitrary direction, at zero temperature. 




hz\\ = VI ~ '^'i ^ 2/3 
h^^ = 2/(33/2a), a > 2/3. 



(20) 



(left). 
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Wc also considered the effect of temperature on tlie Stoner-Wohlfarth model 
but at quasi-equilibrium. Now, we come to deal with the thermoactivatcd 
switching of the particle's magnetic moment, a process that occurs at short- 
time scales. At nonzero temperatures the magnetization vector of the particle 
can surmount the energy barrier due to thermal fluctuations as argued by 
Neel [1]. This effect is particularly pronounced for small particles with lower 
values of the potential barrier (5). Indeed, the magnetization vector of the 
particle shuttles between the two energy minima and the characteristic time 
for this thermoactivated rotation of the spin over the anisotropy barrier A7^ 
is approximately given by the celebrated Van't Hoff-Arrhenius law [14], [15] 

r = Toe^«/'=«'^, (21) 



where Tq is usually taken as temperature and field independent and on the or- 
der of 10"^° — 10~^^s. To is not necessarily the same for different ferromagnetic 
materials. It can indeed be assumed as constant only if the magnetization 
vector is always in one of the energy minima, which happens only if the the 
minima have zero width, or equivalently if the barrier is infinitely high. How- 
ever, in any realistic case, there is a finite probability that the magnetization 
vector spends some time in the vicinity of either minimum, in which case the 
prefactor need not be constant, and certainly depends on temperature and 
field. Indeed, in the case of cubic anisotropy the assumption of a constant 
prefactor tq turns out to be a bad approximation [16]. For a given measur- 
ing time Tm the system is at thermal equilibrium when r„j ^ r, and in this 
case the particle is said to be superparamagnetic, which corresponds to the 
temperature range 

ln(T^/To) > An/kBT > 0. (22) 



In fact, it was argued in [17] that due to the smallness of Tq the range of ther- 
mal equilibrium can extend down to temperatures at which the energy-barrier 
height is much larger than thermal energy. More precisely, for magnetic mea- 
surements with Tm ~ 100s, this range is < AT-C/ksT < 25, which means that 
it is too restrictive to argue that superparamagnetism occurs only when ther- 
mal energy is on the order of the energy-barrier height. In the case r ^ the 
particle's magnetization does not change during the time of observation, and 
the particle exhibits stable ferromagnetism and is said to be in a blocked state. 
The temperature at which such a transition occurs, namely the temperature 
at which the relaxation time r is equal to the observation time r^, is called the 
blocking temperature and is denoted by Tg. For nanoparticle assemblies with 
size distribution, the same temperature may sometimes be above Tb for some 
particles, and below it for the others. Such systems may thus behave as super- 
paramagnetic for some high values of the temperature T, as ferromagnetic at 
low values of T, and as a mixture of both at intermediate T. For such a system 
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a different characteristic temperature is used, this is denoted by Tmax, and is 
roughly given, at least in the dilute case, by an average of all Tg's. One should 
not forget however that the scale of 100s depends on the experimental appa- 
ratus. Indeed, in Mossbauer effect measurements, the experimental time is the 
time of Larmor precession, which is on the order of 10~^s, while in neutron 
scattering experiments it is on the order of 10~^^s. In addition, the time scale 
can be completely different in different areas of applications. For example, in 
magnetic recording, in order to keep the data stored on a magnetic tape for 
years, one should have r ^ lO^s, and in rock magnetism the magnetization 
decays (or relaxes) within geological times which may be millions of years. 

Owing to these important practical applications and many others, the relax- 
ation time of the particle's magnetization is a very important and fundamental 
physical quantity that deserves extensive and rigorous investigation. This ac- 
tually started, in this context, with the work of Kramers on transition-state 
method [18]. Kramers showed, by using the theory of the Brownian motion, 
how the prefactor of the reaction rate (inverse of relaxation time), as a func- 
tion of the damping parameter and the shape of the potential well, could be 
calculated from the underlying probability- density diffusion equation in phase 
space, which for Brownian motion is the Fokker-Planck equation (FPE). He 
obtained, by linearizing the FPE about the potential barrier, explicit results 
for the escape rate for intermediate-to-high (IHD) values of the damping pa- 
rameter and also for very small values of that parameter. Subsequently, a 
number of authors [19], [20] showed how this approach could be extended to 
give formulae for the reaction rate which are valid for all values of the damping 
parameter. These calculations have been reviewed by Hanggi et al. [21]. 

Kramers theory, originally developed for mechanical particles, was first adapted 
to the thermal rotational motion of the magnetic moment by Brown [2] in 
order to improve upon Neel's concept of the superparamagnetic relaxation 
process, which implicitly assumes discrete orientations of the magnetic mo- 
ment and which does not take into account the gyromagnetic nature of the 
system. Brown in his first explicit calculation [2] of the escape rate confined 
himself to axially symmetric (functions of the latitude only) potentials of the 
magneto-crystalline anisotropy so that the calculation of the relaxation rate 
is governed [for potential-barrier height significantly greater than ksT] by the 
smallest non-vanishing eigenvalue Ai of the Sturm-Liouville equation associ- 
ated with the one-dimensional FPE. Thus the rate obtained is valid for all 
values of the damping parameter. As a consequence of this very particular 
result, the analogy with the Kramers theory for mechanical particles only 
becomes fully apparent when an attempt is made to treat non axially sym- 
metric potentials of the magneto-crystalline anisotropy that are functions of 
both the latitude and the longitude. In this context. Brown [2] succeeded in 
giving a formula for the escape rate for magnetic moments of single-domain 
particles, in the IHD damping hmit, which is the analog of the Kramers IHD 
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formula for mechanical particles. In his second 1979 calculation [3] Brown only 
considered this case. Later Klik and Gunther [22], by using the theory of first- 
passage times, obtained a formula for the escape rate which is the analog of the 
Kramers result for very low damping. All these (asymptotic) formulae which 
apply to a general non-axially symmetric potential, were calculated explicitly 
for the case of a uniform magnetic field applied at an arbitrary angle to the 
anisotropy axis of the particle by Geoghegan et al. [23] and compare favorably 
with the exact reaction rate given by the smallest non- vanishing eigenvalue of 
the FPE [24], [25], [26], [27] and with experiments on the relaxation time of 
single-domain particles [26], [28]. 



In both limiting cases considered by Brown the time dependence of the av- 
erage magnetization is a single exponential and the relaxation rate is given 
by the eigenvalue Ai previously mentioned. Subsequently, Ai was numerically 
calculated by Aharoni for arbitrary values of AT^/ZcsT without a magnetic 
field [29] and in a longitudinal magnetic field [30]. The correction terms for 
the high-barrier result for Ai were given by Brown [31]. Various limiting cases 
and the corresponding expressions for Ai were further investigated in [32], 
[33], [34], [35], [24]. Let aside such limiting cases, the FPE for an assembly 
of single-domain ferromagnetic particles cannot be solved analytically. The 
magnetization relaxation curve in fact consists of an infinite number of ex- 
ponentials. In this case, it was shown in Refs. [36], [37], [38] that it is more 
convenient to introduce the so-called integral relaxation time Tint, defined as 
the area under the relaxation curve after a sudden infinitesimal change of the 
magnetic field. Tint depends on all eigenvalues \k,k = 1,2,... and therefore 
contains more information than can be rendered by the first eigenvalue Ai, 
and more importantly, it can be measured experimentally. Moreover, it turns 
out that, unlike Ai, Tint can be calculated analytically for uniaxial particles in 
the longitudinal magnetic field for the arbitrary values of the damping param- 
eter, and Tint recovers the analytical results of Brown for Ai in the asymptotic 
regions. It was shown in [37] that the drastic deviation of t~^1 from Ai at low 
temperatures, starting from some critical value of the field, is a consequence 
of the depletion of the upper potential well. In these conditions, the integral 
relaxation time T^t consists of two competing contributions corresponding to 
the overbarrier and intrawell relaxation processes. 



After the above survey, we now give some of the most commonly used ex- 
pressions for the relaxation time, together with some cautionary remarks. As 
indicated above. Brown [2] at first derived a formula for Ai, for an arbitrary 
axially symmetric bistable potential of Fig. 1 with minima at = (0, vr) sepa- 
rated by a maximum at 9m = arccos(— /i), which when applied to Eq. (3) for 
hjje, i.e. a magnetic field parallel to the easy axis, leads to the form given by 
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Aharoni [30], 

Ai ~ _ h') X [(1 + h) e-'^^(^+'^)' + {l-h) e-^^'^'-'^n , (23) 



where 7 is the gyromagnetic factor, a the damping parameter. < h < 1, 
h — 1 being the critical value at which the bistable nature of the potential 
disappears [see section 2.1]. 

In order to describe the non-axially symmetric asymptotic behavior, let us 
denote by /3A7^_ the smaller reduced barrier height of the two constituting 
escape from the left or the right of a bistable potential [see Eq. (5) for the 
axially-symmetric case]. Then for very low damping, i.e. for a x /5A7i_ ^ 1 
(with of course the reduced barrier height pATi- ^ 1, depending on the size 
of the nanoparticle) , we have [2], [27] the following asymptotic expression for 
the inverse relaxation time 



TvId ^ ^ (24) 
c^^{ujix PiHo - 7ii)e-^^^°-^i) +UJ2X PiHo - 7^2)6-^^^°"^^)} , 
where a is the damping parameter and 

a 27 



is the free-diffusion time, i.e., the characteristic time of diffusion in the ab- 
sence of the potential. For the IHD limit, where a x /3A7Y_ > 1 (again with 
the reduced barrier height /5A7i_ much greater than unity) we have [27] the 
asymptotic expression 

rf^^ ^ {a;,e-^(^°-^^) + u^e-^^-^^-^^^} , (26) 



where 



m 



(i).(i) 

2 5 



2^1 



UJ2 



m 



(2)^2) 
2 5 



2^1 



m 



(0)^(0) 

2 ) 



2C1 'c. 



2m 



4cS°>cf 



Here uji^uj2 and uoq are respectively the wells and saddle angular frequencies 
associated with the bistable potential, VLq is the damped saddle angular fre- 
quency and the cf' are the coefficients of the truncated (at the second order 



16 



in the direction cosines) Taylor series expansion of the crystalhnc anisotropy 
and external field potential at the wells of the bistable potential denoted by 1 
and 2 and at the saddle point denoted by 0. A full discussion of the applica- 
tion of these general formulae to the particular potential, which involves the 
numerical solution of a quartic equation in order to determine the cf^ with 
the exception of the particular field angle ^0 = | or |, in Eq.(l) is given in 
Refs. [23], [39]. 

In Ref. [40], we used either Eqs. (23) and (24) or (23) and (26), and solved 
the equation r — Tm for the blocking temperature Tg, as a function of the 
applied field, for an arbitrary angle between the easy axis and the applied 
magnetic field. In particular, for very small values of %p we used Eq. (23), as the 
problem then becomes almost axially symmetric and the arguments leading to 
Eqs. (24) and (26) fail [2], [3], [23], and appropriate connection formulae had 
to be used so that they may attain the axially symmetric limit. Application of 
the above asymptotic formulae to the calculation of the blocking temperature 
Tb (or Tmax) as a function of the applied field appeared to recover the experi- 
mental observation, but this result turned out to be spurious. An explanation 
of this behavior follows (see also [41], [42]): in the non-axially symmetric IHD 
asymptote (26), which is formulated in terms of the Kramers escape rate, as 
the field tends to zero, for high damping, the saddle angular frequency ujq 
tends to zero. Thus the saddle becomes infinitely wide and so the escape rate 
predicted by Eq. (26) diverges [see similar discussion in Sect. 3] leading to 
an apparent rise in the blocking temperature until the field reaches a value 
sufficiently high to allow the exponential in the Arrhenius terms to take over. 
When this occurs the blocking temperature decreases again in accordance 
with the expected behavior. This is the field range where one would expect 
the non-axially asymptote to work well. In reality, as demonstrated by exact 
numerical calculations of the smallest non vanishing eigenvalue of the Fokker- 
Planck matrix, the small field behavior is not as predicted by the asymptotic 
behavior of Eq. (26) (it is rather given by the axially-symmetric asymptote) 
because the saddle is limited in size to loq. Thus the true escape rate cannot 
diverge, and the apparent discontinuity between the axially-symmetric and 
non axially-symmetric results is spurious, leading to an apparent maximum in 
Tsih). In reality, the prefactor in Eq. (26) can never overcome the exponential 
decrease embodied in the Arrhenius factor. Garanin [41] has discovered bridg- 
ing formulae which provide continuity between the axially-symmetric Eq. (23) 
and non-axially symmetric asymptotes leading to a monotonic decrease of the 
blocking temperature with the field in accordance with the numerical calcu- 
lations of the lowest eigenvalue of the Fokker- Planck equation and also in 
agreement with the fact that the field suppresses the energy barrier upon 
which the blocking temperature decreases. 

An illustration of this was given in Ref. [41] [see also [42] for more detail] for 
the particular case oi t/j — 7r/2, that is a transverse field. If the escape rate is 
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written in the form 

= -Aexp{-f3An) 

TT 

where k is the attempt frequency given by 
m 

then the factor A, as predicted by the IHD formula, behaves as a/Vh for 

h <S 1,q;^, while for h — 0, A behaves as 2TrayJa/Tr, which is obviously dis- 
continuous. So, a suitable interpolation formula is required. Such a formula 
(analogous to that used in the WKBJ method [43]) is obtained by multiplying 
the factor A of the axially symmetric result by e~^/o(0; where /o(0 is the 
modified Bessel function of the first kind, and ^ = 2ah [see [41]]. This interpo- 
lation formula, as is obvious from the large and small ^ limits, automatically 
removes the undesirable 1 / \/h divergence of the IHD formula and establishes 
continuity between the axially symmetric and non-axially symmetric asymp- 
totes for ip = 7t/2, as dictated by the exact solution. 

It is apparent from the discussion of this section that the Neel-Brown model 
for a single particle is unable to explain the maximum in Tmax, observed in 
experiments, as a careful calculation of the asymptotes demonstrates that they 
always predict a monotonic decrease in the blocking temperature. A possible 
explanation of the maximum of T^ax as a function of the applied field was 
given in [40] where it was shown that the non-linearity of the field dependence 
of the superparamagnetic contribution to the assembly magnetization, volume 
distribution, and anisotropy are responsible for this effect. 



3 A first step beyond the Neel-Brown model: effect of exchange 
interaction 

As discussed in the introduction, in order to calculate the relaxation time 
of a nanoparticle in an inhomogeneous magnetic state induced by surface ef- 
fects, one has to resort to microscopic theories. Then, one has to take account 
explicitly of microscopic interactions, such as exchange and/or dipolar inter- 
actions, in addition to the magneto-crystalline and surface anisotropy and 
applied magnetic field. However, the problem associated with the generaliza- 
tion of Brown's theory to include interactions and eventually surface effects 
is really involved and can in general only be tackled numerically. But before 
one can attack this problem, one needs to understand the effect of, e.g., ex- 
change interaction on the relaxation time of the minimal system, i.e., a pair of 
atomic spins coupled via exchange interaction, including of course the usual 
magneto-crystalline anisotropy and Zeeman terms. Besides, this is the unique 
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non-trivial step towards the above-mentioned generalization, where analytical 
expressions can be obtained for the relaxation time. It was the purpose of 
Ref. [44] to solve this problem and to compare with the Neel-Brown result 
in the one-spin approximation. Accordingly, we studied the effect of exchange 
coupling on the relaxation rate of a magnetic system of two spins within 
hanger's approach. We found a particular value of the exchange coupling, 
that is j = J / K = jc = 1 — h"^, where h = H/2K, at which the number of 
saddle points changes. For j < there are two saddle points and the reversal 
of the spins proceeds in two steps: the first spin crosses one of the saddle points 
into the anti-ferromagnetic state, and then the other spin follows across the 
second saddle point, to end up in the ferromagnetic order of opposite direc- 
tion with respect to the initial one. For j > jc, the two spins are so strongly 
exchange-coupled that they cross a single saddle point together. For j ^ 
the Neel-Brown result for a single spin is recovered. As a byproduct, we also 
showed that Langer's quadratic approximation of the potential energy at the 
saddle point fails when the exchange coupling assumes the critical value jc 
even in the IHD limit. 

It is worth mentioning that the effect of dipolar interactions on the relaxation 
time of magnetic systems has been studied by a few authors. In Ref. [45], 
[46] this was done for assemblies of nanoparticles where approximate expres- 
sions were obtained using the (mean-field) barrier-height approach that is only 
valid in the limit of infinite damping (no gyromagnetic effect). In Ref. [47] the 
authors used perturbation theory and obtained an expression for the relax- 
ation rate in the case of weakly interacting superparamagnets with various 
anisotropy orientations and in the absence of external field. Finally, the au- 
thors of Ref. [48] dealt with a pair of coupled dipoles using the (numerical) 
Langevin approach and found some results similar to ours [see below]. 

In this section we will first give a brief account of Langer's approach [49] (see 
also [51] for uniform (one-spin problem) and non-uniform magnetization and 
[52] for comparison with Kramers' theory). 



3.1 Switching rate in Langer's approach 

Within this approach the problem of calculating the relaxation time for a 
multidimensional process is reduced to solving a steady-state equation in the 
immediate neighborhood of the saddle point that the system crosses as it goes 
from a metastable state to another state of greater stability. The basic idea 

[49] is that "... one imagines setting up a steady-state situation by continuously 
replenishing the metastable state at a rate equal to the rate at which it is leaking 
across the activation- energy barrier. By identifying the current flowing over 
the barrier with the desired decay rate, one avoids having to solve the complete 
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time- dependent problem..." , especially in the multidimensional case where this 
problem is too difficult to solve. However, this approach is only valid in the 
limit of intermediate-to-high damping because of the assumption, inherent to 
this approach, that the potential energy in the vicinity of the saddle point may 
be approximated by its second-order Taylor expansion (see below) . The result 
for small damping fails because the region of deviation from the Maxwell- 
Boltzmann distribution, set up in the well, extends far beyond the narrow 
region at the top of the barrier in which the potential may be replaced by its 
quadratic approximation. Therefore, in Langer's approach, one concentrates 
on the distribution function p{{'r]},t) as the probability that the system is 
found in the configuration {77} at time t. The time evolution of p is governed 
by the following equation 

dp^fdp\ ^fdp\ ^27) 

dt y dt J \ dt J j^ly^f,^ 

where the first term accounts for the "internal" dynamics of the system de- 
scribed by 

^ = EA.|^, (28) 
dt Y % 



where A is a, fully anti-symmetric matrix. For a magnetic system, as in our 
case, the analog of Eq. (28) is given by the Landau-Lifshitz equations (50). 
The second term in (27) accounts for the dynamics of the system driven by 
its interaction with the heat bath, and is given by 



where j3 = l/ksT and the constant Fj characterizes the variation rate of f]i{t). 
Combining Eqs. (28) and (29) in Eq. (27) leads to the Fokker-Planck equation 

dtp + Yldr„Ji-Q, (30) 

i 

which is a continuity equation with the probability current 



j,^-Y:M,,(d,^H+^d,^^p, 



(31) 



and where the matrix M reads 

Mij = f^ViSij-Aij. (32) 
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In order to calculate the nucleation rate, one must solve Eq. (30). A particular 
solution is obtained at equilibrium and is given by the Maxwell-Boltzmann 
distribution 

Peg = exp(-/37i{?7})/Zo, (33) 



which corresponds to zero current, Jj = 0. However, what is really needed 
is to obtain a finite probability current fiowing across the saddle point r]^^\ 
As stated above, considering a steady-state so that the rate at which the 
metastable state is replenished is equal to the rate of leak across the activation- 
energy barrier, the problem then consists in finding this steady-state as a 
solution of Eq. (30) in the immediate neighborhood of the saddle point. To 
obtain this solution, it is more convenient to work in the frame of canonical 
coordinates with their origin at rj^^\ defined by 

^n^ET^ni{Vi-vh: (34) 



where the transformation matrix V diagonalizes the energy Hessian, so that 
the energy can be rewritten near {^} = {0} as 



1 TV 

^ ) + - j: Kil (35) 



n=l 



The A„'s are the eigenvalues of the energy Hessian. By the definition of ri^^\ 
one of the A's, denoted in the sequel by A_, is negative; that is, the energy at 
the saddle point diminishes on either side of the surface ^_ = 0. Moreover, we 
expect that m < N eigenvalues vanish, which corresponds to the fact that the 
saddle point does not possess all of the symmetries of 7i {r/}. The unbroken 
symmetries produce Goldstone modes ^ that must be eliminated from the 
nucleation rate as they cause the latter to diverge. Consequently, the saddle 
point is actually a bounded, finite-dimensional subspace of the ?7-space. The 
volume of this subspace will be denoted below by V. For the single-spin prob- 
lem, this is just the space spanned by rotations around the angle ip, while 
cos 6"-*^ = — /i, and its volume is given by = 27rsin^^*^*^ = 27rVl — h'^- 
Next, one diagonalizes the transition matrix that results from the dynami- 
cal equation (28) linearized at the saddle point, that is one has to solve the 
eigenvalue problem: 

K M^m Um = I^Un, (36) 

m 

^ A spontaneous breaking of a continuous symmetry entails the existence of a mass- 
less mode, that is a zero-energy fluctuation, called the Goldstone mode: This is the 
well-known Goldstone theorem [50]. 
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where the prime on the summation indicates that we omit all n's for which 
A„ = 0, and 

= V^ii-Aij)V^j = -M^rn- (37) 



As was argued by Langer [49], one, and only one, of the eigenvalues k, must be 
negative. Indeed, if the saddle point is to describe the nucleating fluctuation, 
there must be exactly one direction of motion away from the saddle point 
{^} = {0} which the solution 

(U^^ne-'^*, (38) 
of the equation of motion of the ^„ modes 



dt 



J^^nrnXmiU): (39) 



is unstable. 

Finally, Zq in Eq. (33) is a normalization factor defined by the condition that 
the equilibrium probability density (33) is normalized at the metastable state, 
i.e., 




Vr] exp(-/??i("^){r7}) = 1, 



which leads to 




where A™ are the eigenvalues of the energy Hessian at the metastable state 
and E^^^ its energy. 

Taking all this into consideration, Langer has given a formula for the switching 
rate T out of the metastable state which reads. 




where AEq = Eq — Eq is the barrier energy between the saddle point and 
the metastable state, and the double prime on the product symbol indicates 
that we omit the negative eigenvalue A_ and all n for which A„ = 0. 
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The whole procedure can be summarized as follows: From a static study of the 
energy, one identifies the metastable states and saddle points. Then at a given 
saddle point one expands the energy to second order, calculates the Hessian 
and its eigenvalues. The vanishing eigenvalues corresponding to the Gold- 
stone modes, associated with the unbroken symmetries at the saddle point, 
are eliminated from the final expression of the relaxation rate. The volume V 
of the saddle-point subspace is also calculated. Then, one solves the eigenvalue 
problem (36) and obtains the negative eigenvalue k. Finally, one calculates the 
eigenvalues of the Hessian of energy expanded at the metastable state. 

A quite similar approach has been used by D.A. Garanin [53] to calculate 
the thermo-activation rate of charged particles in a magnetic field. Here, the 
Fokker-Planck equation is written in terms of the 6 variables, 3 spatial coordi- 
nates and 3 momenta, and then linearized near the saddle point. The problem 
then amounts to computing the flux and normalization in the wells rendering 
a compact expression for the rate in terms of the frequencies in the wells and 
at the saddle point. 



3.2 Relaxation rate for the two-spin problem within hanger's approach 

We consider a system of two (classical) exchange-coupled spins with the Hamil- 
tonian 



where J > is the exchange coupling, K > Q the anisotropy constant, H the 
applied magnetic field, and h the reduced field, i.e., < < 1. Cj are uniaxial 
anisotropy unit vectors. We use the dimensionless units, i.e., \H] = 2K for 
energy and [t] = 1/(27 A') for time. Here we restrict ourselves to the case 
h II ei,i = 1,2, with the same anisotropy constant. Owing to the symmetry 
of this system with respect to rotations around the easy axis, the number of 
variables reduces to three, 9i,92,(fi = (pi — (fi2- 

Now we apply hanger's approach to the energy (42) and study the relaxation 
rate as a function of the exchange coupling j . We first analyze the energyscape 
in Fig. 5. The absolute minimum of the energy (42) corresponds to the ferro- 
magnetic ordering of the two spins along the easy axis. 



(si ■ ei)^ + (s2 • 62)^ -h - {si + S2), 



j = J/K, h = H/2K 



(42) 



(0,0, (^); e(°) = -|- 2/1-1, 



(43) 
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Fig. 5. Energyscape from Eq. (42) {(p = 0) for j = 0.25,0.8, 1.5 and h = 0.1. The 
arrows indicate the switching paths. 



where henceforth e'^'^^ denotes the energy of the state. One metastable state 
corresponds to a ferromagnetic ordering opposite to the field 

(7r,7r;^); e^^\ = + 2h - 1. (44) 

There is also the metastable state of anti-ferromagnetic ordering, 

(0,7r;(^)or(7r,0;(^); e^^) = | - 1. (45) 

As to saddle points, we find that their number and loci crucially depend on 
the exchange coupling constant j. More precisely, for j > jc = 1 — h'^ , there is 
a single saddle point given by 

(cos^i = COS ^2 = -h;(py, ef^ = + ^^ (46) 
whereas for j < jc there are two saddle points given by 

cos ,2 = \(^-h-a'± VA?) = X^; ^, (47) 
where £ = ± and 
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A' ^ {h + a'f + 2j -4b'h, 



urn' 



1-i 



with energy, 




(48) 



At j — jc the saddle point (47) with e — + merges with the saddle point 
(46), while that with e — — merges with the metastable state (45), see Fig. 5 
central panel. 

Starting with both spins aligned in the metastable state (44) with 9i^2 = tt, if 
j < jc one of the two spins crosses the saddle point (47) {e — +) into the state 
(45) by reversing its direction. Then the second spin follows through the second 
saddle point (47) (e = — ), of lower energy due to the exchange coupling (see 
Fig. 6), ending up in the stable state (43). In Fig. 5 (left) the path is indicated 
by a pair of curved arrows. There are actually two such paths corresponding 
to the two-fold symmetry of the problem owing to the full identity of the two 
spins. Note that when the first spin starts to switch and arrives at 9i ~ 7r/2, 
the second spin has ^2 ^ tt (hence the curved arrows in Fig. 5), which suggests 
that in the switching process of the first spin, the position of the second spin 
undergoes some fluctuations creating a small transverse fleld, and when ^1 = 
the second spin goes back to the position 62 — n before it proceeds to switch 
in turn. The successive switching of the two spins through the corresponding 
saddle points is a sequential two-step process ^ , so the relaxation rates for 
j < jc add up inverse- wise. In the case j > jc the two spins cross the unique 
saddle point (46) to go from the metastable state (44) into the stable one 
(43) in a single step, see Fig. 5 (right) where the path is indicated by a single 
straight arrow. There is the symmetry ±d^^\ which leads to a factor of 2 
in the relaxation rate. Therefore, if we denote by ^^<j^, ^J<jci ^j>jc 
respective relaxation rates, the relaxation rate of the two-spin system can then 
be written as 



^ Similarly, it was found in Ref. [48] that the reversal of the two dipoles considered 
is a two-stage process with an intermediary metastable antiparallel state. 



^J<jSj<J(^J<Jc+'^J<jJ^^j<jc 



r = <^ 



(49) 
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j>jc1 



if j > jc- 
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Fig. 6. Energy barrier as a function of j. 

In Langer's approach the F's in Eq. (49) arc obtained from the steady-state 
Fokker-Planck equation (SSFFE) hnearized around each saddle point, using 
the fluctuating variables rji = {ti,Pi) where ti = 9i — 9^, = (pi — ipf, i = 1,2, 
or more adequately the "canonical" variables ■0n — (C±, C±) with ^± = {ti ± 
t2)/V^: C± — (Pi ='=P2)/'\/2, in which the energy Hessian is diagonal. 

The deterministic dynamic of the system is governed by the Landau-Lifshitz 
equations, which upon linearization near the saddle point, read 

dtU — —d^Ho — adt-Ho 

< * * P. 2 t, 2 ^^^^ 



where 7^2 is the quadratic approximation of the energy (42) near the saddle 
point and a is the damping parameter. In matrix form Eqs. (50) become 

dtVi^J2Mijdr,,n2, (51) 



where M is the dynamic matrix containing the precessional and dissipative 
parts. Rather than investigating the stochastic trajectories r)i{t) that arise by 
adding a noise term in Eq. (51), Langer concentrates on the distribution func- 
tion p{{ri},t) as the probability that the system is found in the configuration 
{rj} at time t. The time evolution of p is governed by the Fokker-Flanck equa- 
tion (30), where now the current Jj is given by Eq. (31) with the Hamiltonian 
Ti. replaced by its quadratic approximation H2- 
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Instead of solving the time-dependent equation (30), Langer solves the SSFPE 
dtp = near the saddle point. The steady-state is realized by imposing the 
boundary conditions: p ~ Peg near the metastable state and p ~ beyond 
the saddle point. Then, the problem of calculating the escape rate boils down 
to the calculation of the total current by integrating the probability current 
(31) over a surface through the saddle point. After performing all these steps, 
Langer arrives at his famous expression (41) for the relaxation rate which is 
valid in the IHD limit. The general calculations can be found in greater details 
in [49], or applied to magnetic particles in [51], see also the review article [52]. 

However, since the escape rate is simply given by the ratio of the total current 
through the saddle point to the number of particles in the metastable state, it 
turns out that hanger's result for the escape rate can be retrieved by only com- 
puting the energy-Hessian eigenvalues near the saddle points and metastable 
states, from which one then infers the partition function Zs of the system 
restricted to the region around the saddle point where the energy-Hessian 
negative eigenvalue is (formally) ^ taken with absolute value, and the parti- 
tion function of the region around the metastable state. When computing 
these partition functions, one has to identify and take care of the Goldstone 
modes associated with the unbroken global symmetries of the different states. 
Finally, one computes the unique negative eigenvalue k, of the SSFPE corre- 
sponding to the unstable mode at the saddle point. More precisely, k, is given 
by the negative eigenvalue of the dynamic matrix M„i„ = —Xn{VMV'^)mn, 
where the A„'s are the eigenvalues of the Hessian at the saddle point and V 
is the transformation matrix from rji to ipn [see Eq. (34)]. 

Consequently, Langer's final expression for the escape rate is rewritten in the 
following somewhat more practical form 



where \k\ is the attempt frequency which contains the damping parameter a. 



3.3 Relaxation rate of the two-spin system 

Now, we give the different relaxation rates for j > jc, j < jc, and j ~ jc- 



^ See Ref. [49] for a rigorous derivation 
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3.3.1 Relaxation rate for j > jc 



The quadratic expansion of tlie energy (42) at tiie saddle point (46) reads, 

HP = - ie + + f e, (53) 



with zero eigenvalue for the C+ mode, that is the Goldstone mode associated 
with the rotation around the easy axis, e^*^^ is given in Eq. (46). The negative 
eigenvalue of the SSFPE, corresponding to the unstable mode, is k, — —ajc- 

The partition function around the saddle point Zg = J dVlidVL2e~^^'' , where 
dVLi = sinOidOidipi, i = 1,2 is calculated by changing to the variables C±,C±, 
setting sin 9i ~ ^/J^ (see Eq. (46)) in the integration measure, and substituting 
27r for the integral over and finally computing the Gaussian integrals. 
Hence, 




where we have formally replaced the negative eigenvalue {—jc) by its absolute 
value. Next, the partition function at the (well defined) mctastable state 
(44) is computed by expanding the energy up to 2"*^ order, leading to 

Z ^e-^^'^U-X ^ (55) 

[pj ii-h)ij + i-hy ^''^ 



Finally, using Eq. (52), the relaxation rate for j > jc reads, upon inserting the 
symmetry factor of 2, 



2jA(l_/,2)(i_^)l + (l 

^l-jc/j 



Ae(°) = (1 - h)^. 



(56) 



It is readily seen that for j — > oo (56) tends to the Neel-Brown result. 



r,>,,^2jA(l-/.^)(l-/.) e-^^-'^)^ 



(57) 



for the relaxation rate of one rigid pair of spins with a barrier height twice 
that of one spin. In fact, the result (57) coincides with the Neel-Brown result 
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in Eq. (23) upon reinstating the energy and time units and remembering that 
the result (57) is for the one-way barrier crossing. Note that the convergence 
of (56) to (57) is so slow that Tj^j^ remains above the Neel-Brown result 
and only merges with it for j > 10. This confirms the fact that the one-spin 
approximation is only valid for extremely large exchange interaction. 



3.3.2 Relaxation rate for j < jc 

Here the attempt frequencies |«;^| (for e — ±) cannot be obtained in a closed 
form and are thus computed numerically, this means that the present case can 
only be dealt with semi-analytically. On the other hand, following the same 
procedure as for j > jc we obtain the relaxation rate for £ = ±, 




where Ae^°^ is the barrier height given by the energy (48) measured with 
respect to (44) or (45) for e — +,— respectively, and (see Eq. (47) et seq. for 
notation) 



= ^{1 + Ifh - j/2Y - (a- + h)\ 

N+^{l-h){j + l-h), N-^jc-j. 

The hmit of the relaxation rate (58) when j — > is just the Neel-Brown 
result for one spin. Indeed, the product of the last two factors in the prefactor 
tend to (1 — h)/\/2, the attempt frequencies tend to ajc — a{l — h?), and the 
energy barriers Ae(°) ^ {l-hf/2. 

Note that in the present regime of j < jc the large value of anisotropy has not 
changed the temperature dependence of the individual relaxation rates, i.e., 
^)<jc^ with e = ±. This is due to the fact that 1/ \/T appears in the prefactor 
each time there is a continuously degenerate class of saddle points [51], which 
is indeed the case for j > jc and j < jc with e — + and e — —. However, 
anisotropy do affect the temperature dependence of the relaxation rate of the 
two-spin system, since for j < jc there are two saddle points bringing each a 
factor 1/\/T, see the first line in Eq. (49). 



3.3.3 The case of j ~ jc 

When j approaches jc either from above or from below, more Hessian eigen- 
values (in addition to A^+) vanish, rendering the saddle point rather fiat and 
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thus leading to a divergent relaxation rate. Indeed, for j ~ jc the relaxation 
rate (56) diverges, which clearly shows that hanger's approach making use of 
a quadratic approximation for the energy at the saddle point, e.g., Eq. (53), 
fails in this case. The remedy is to push the energy expansion to the 6*^— order 
in the variable ^_ (since Ag_ = j — jc — ^ as j — > jc), i.e., 

Ses = e, - ei°) ^ + + ^C, (59) 



where 

C = — < 0, = > 0. 

3 30 



Then, the contribution y2'7r//5Ag_ of the mode ^_ to the relaxation rate (56) 
must be replaced hjJ^^dC,^e~'^^^, upon which the divergence of Tj^j^ is cut 
out (see Fig. 7). 




Fig. 7. Relaxation rate (in logarithmic scale) of the two-spin system in the case 
of intermediate-to-high damping. The horizontal line is the Neel-Brown result with 
the barrier height taken twice that of a single spin. 

Similarly, when j approaches jc from below, for e = — the eigenvalue A^- 
vanishes at j = jc, upon which the saddle point (47) {s = —) merges with 
the state (45) and thereby the partition functions Zs and Zm tend to infinity. 
However, as Zm oc l/(jc — j) increases much faster than Zg, and 0, 
^J<jc tends to zero as j — > jc- On the other hand, for e = +, both A^- and A^- 
vanish leading to a divergent relaxation rate, since now Zg diverges but Zm in 



30 



Eq. (55) remains finite for tlie metastable state (44) is well defined. Indeed, 
as j — > jc, the relaxation rate r^<j^ goes over to the result in Eq. (56) upon 
making the change j jc, and taking account of the symmetry factor. In this 
case, the divergence at the point j = jc cannot be cut ofT by expanding the 
energy beyond the 2"*^ order and r^^^^ is simply cut off at the point where it 
joins ^j>j^ taking account of Eq. (59). 

In Fig. 7 we plot {In of) the relaxation rate of the two-spin system as defined 
in Eq. (49), under the condition of IHD, in which Langer's approach is vahd, 
that the reduced barrier height /^Ae^^^ » 1 and af3Ae^^^ > 1 [52]. We also 
plot separately both relaxation rates for j < jc- We see that the relaxation 
rate of the two-spin system contains two unconnected branches corresponding 
to the two regimes, j < jc and j > jc, the bridging of which would require 
a more sophisticated approach. Fig. 7 also shows that as j increases, but 
3 <S jc, the relaxation rate r^<^-^ decreases because the switching of the first 
spin is hindered by the (ferromagnetic) exchange coupling. While ^J<j^ is an 
increasing function of j with a faster rate, since now the exchange coupling 
works in favor of the switching of the second spin. This is also illustrated by the 
evolution of the energy barrier height in Fig. 6. As j approaches jc from below, 
the relaxation rate ^'j<j^ tends to Tj>j^ because the respective saddle points 
merge at j = jc- Whereas ^J<j^ goes to zero since the corresponding saddle 
point merges with the antiferromagnetic state that is no longer accessible to 
the system. For j > jc, as j increases the minimum (43), the metastable state 
(44) and the saddle point (46) merge together, which means that the system is 
found in an "energy groove" along the direction 9i — 02 because the eigenvalue 
corresponding to the mode 61 — 62 becomes very large, and thus the escape 
rate decreases and eventually reaches the Neel-Brown value at large j. 



4 Discussion: What is wrong with the one-spin approximation ? 

We have seen that the magnetization of a nanoparticle can overcome the en- 
ergy barrier and thus reverse its direction, at least in two ways: either under 
applied magnetic field which suppresses the barrier, or through thermally ac- 
tivated statistical fiuctuations. Within the framework of the one-spin approxi- 
mation, the magnetization switching under applied field, at zero temperature, 
is well described by the Stoner-Wohlfarth model. At finite temperature and 
short-time scales, crossing of the energy barrier activated by thermal fiuc- 
tuations is described by the Neel-Brown model and its extensions reviewed 
above. Both of these models have been confirmed by experiments on individ- 
ual cobalt particles [6]. At finite temperature, but at quasi-equilibrium, the 
magnetization switching occurs according to two distinct regimes. At very low 
temperature, this is due to the coherent rotation of all spins, as in the Stoner- 
Wohlfarth model, whereas at higher temperature, the magnetization switches 
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by changing its magnitude. This results in a shrinking of the Stoner-Wohlfarth 
astroid as described by the modified Landau theory [7] , and confirmed by ex- 
periments [6]. However, it is clear that the change of magnetization magnitude 
cannot be explained in the framework of the one-spin approximation. Indeed, 
it can only be understood as the result of a successive switching of individ- 
ual (or clusters of) spins inside the particle, which is necessarily a multi-spin 
system. Indeed, deviations from the single-spin approximation, and thereby 
from both the Stoner-Wohlfarth and Neel-Brown models, have been observed 
in metaUic particles [54], [55], and ferrite particles [56], [57]. These deviations 
have materialized in terms of the absence of magnetization saturation at high 
fields, shifted hysteresis loops after cooling in field, and enhancement at low 
temperature of the magnetization as a function of applied field. The latter ef- 
fect has been clearly identified in dilute assemblies of maghemite particles [58] 
of 4 nm in diameter. In addition, aging effects have been observed in cobalt 
single particles and have been attributed to the oxidation of the sample sur- 
face into antiferromagnetic CoO or NiO [see Ref. [6] and references therein for 
a discussion of this issue]. It was argued that the magnetization reversal of a 
ferromagnetic particle with antiferromagnetic shell is governed by two mech- 
anisms that are supposed to be due to the spin frustration at the core-shell 
interface of the particle. On the other hand, according to Mossbauer spec- 
troscopy some of the above-mentioned novel features are most likely due to 
magnetic disorder at the surface which induces a canting of spins, or in other 
words, an inhomogeneous magnetic state inside the particle. 

As argued earlier, understanding these effects requires the development of 
microscopic theories capable of distinguishing and accounting for the various 
crystallographic local environments that develop inside a nanoparticle and 
on its surface. Elaboration of such theories is faced with tough non-linear 
A^-body problems whose study can only be, in principle, efficiently performed 
with the help of numerical approaches. It is however desirable, and indeed even 
necessary, that the numerical calculations be backed by analytical expressions, 
in some limiting cases at least. 

Accordingly, the study of section 3 has helped us understand the effect of 
exchange coupling on the relaxation rate of the two-spin system. The corre- 
sponding analytical results will be very helpful in a generalization to multi-spin 
small particles at least for small values of surface anisotropy and thereby small 
deviations from coUinearity, where it has been shown [59] that the surface 
contribution to the macroscopic energy has a simple cubic anisotropy. How- 
ever, this generalization can only be performed using numerical techniques. 
This is now attempted by the help of the ridge method [60] for probing the 
potential energy surface and locating the saddle points, and by the (Onsager- 
Machlup) path integrals [61] for determining the most probable paths con- 
necting a metastable state to a more stable state. 
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